Non-Markovian large amplitude motion and nuclear fission 

V. M. Kolomietz and S. V. Radionov 

Institute for Nuclear Research, 03680 Kiev, Ukraine 
(Dated: December 8, 2008) 

Abstract 

The general problem of dissipation in macroscopic large-amplitude collective motion and its 
relation to energy diffusion of intrinsic degrees of freedom of a nucleus is studied. By applying 
the cranking approach to the nuclear many body system, a set of coupled dynamical equations for 
the collective classical variables and the quantum mechanical occupancies of the intrinsic nuclear 
states is derived. Different dynamical regimes of the intrinsic nuclear motion and its consequences 
on time properties of collective dissipation are discussed. The approach is applied to the descant 
of the nucleus from the fission barrier. 
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I. INTRODUCTION 



Nuclear large scale dynamics (nuclear fission, heavy ion collisions etc.) is a good probe 
for the investigation of complex time evolution of finite Fermi systems. The conceptual 
question is here how collective modes of motion appear in a system with many degrees of 
freedom and how they interact with all other intrinsic modes. Nuclear collective dynamics 
can be studied by using the concept of macroscopic motion for a few collective degrees of 
freedom, which are chosen to describe gross properties of the nucleus l|, 2], Q]. Such a 
kind of approach is acceptable for a slow collective motion where the fast intrinsic degrees 
of freedom exert forces on the collective variables leading to a transport equation. The 
crucial point of such an approach is the separation of the total energy of the system into 
potential energy, collective kinetic energy and dissipation energy. Moreover the dissipation 
of collective motion implies fluctuations in the corresponding collective variables, as follows 
from the fluctuation-dissipation theorem ^. 

Dissipation of the nuclear collective energy reveals itself, for instance, as the non-zero 
contribution of the internucleonic collisions to the widths of the nuclear giant multipole 
resonances. On the other hand, the experimental observation of the finite variance of the 
kinetic energy of the fission fragments manifests the fact that fluctuations have to be also 
associated with the collective variables. Both the dissipation and the fluctuations can be 
described by the introduction of friction and random forces, related to each other by the 
fluctuation-dissipation theorem. In this respect, the Fokker-Planck or Langevin approaches 
can be used to study the nuclear large scale dynamics, see Refs. [5,6] and references therein. 
In general, basic equations of motion for the macroscopic parameters, describing complex 
dynamics of the many-body systems like nuclei, have non-Markovian structure. One of the 
first considerations of memory (non-Markovian) effects for classical liquids can be found in 
Ref. I?!. For the dynamics of nuclei, memory effects have been investigated within different 



approaches. In this respect, one can mention the d issip ative diabatic model [8!], the linear 
response theory {sj and the fluid dynamic approach [lofQ. In this paper, we would like to 
apply the non-Markovian dynamics to the study of the nuclear fission's characteristics and 
clarify the role of the fluctuation and memory effects in the descant of the nucleus from the 
fission barrier to the scission point. 

The plan of the paper is as follows. In Sect. II we derive the non-Markovian Langevin 
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equation of motion for the nuclear shape variables starting from the collisional kinetic equa- 
tion. Sect. Ill is devoted to details of the numerical determination of the saddle-to-scission 
time and pre-scission kinetic energy in presence of the memory effects and the random force 
for the descent of the nucleus from the barrier to the scission point. Summary and conclu- 
sions are given in Sect. IV. We assume that dynamics of nuclear many body system can be 
described as a coupled motion of macroscopic collective modes and intrinsic nucleonic ones. 
The slow collective modes of the nuclear large-amplitude motion are treated in terms of a set 
of classical time-dependent variables q(t) = {qi{t),q2{t), ...,q]\r(t)}, while the fast intrinsic 
modes are described quantum mechanically through the time-evolution of occupancies of 
nucleonic many body states. 

The intrinsic dynamics can be determined through the Liouville equation for the density 
matrix operator p, 

^ + ^Lmt) = 0, (1) 

where L is the Liouville operator defined in terms of the commutator, Lp = H, p /h of the 
nuclear many body Hamiltonian H{q[t\). 



Using Zwanzig's projection technique and a basis of adiabatic eigenenergies and 
eigenf unctions of the nuclear many body Hamiltonian H{q[t]) [l^, we can get dynamical 
equations for a non-diagonal part of the density matrix, 

Pnm{t) = / dt'qi{t') ^^^^ ^^nm{t tj] [/^ .^^^ (^')p^^ _ (t')p„„ (t')] , (2) 

and its diagonal part, 

X [Pmmit') - Pnnit')] • (3) 

Here, Unm = {En — Em)/fi and matrix elements hi^nm = \dH / dqi\^^ measure the coupling 
between the quantum nucleonic and the macroscopic collective subsystems. 
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To proceed further, we use a random matrix theory developed in our previous paper 
for the case of a single collective coordinate. We omit all intermediate steps and give a basic 
diffusion-like equation of motion for the ensmble averaged occupancies p{E, t) of the many 
body states E = E^- 



where Sij and are, correspondingly, the strengths and widths of the energy distributions 
of the ensemble averaged coupling matrix elements /ii,„m, {(l{t) , <l{t') are the correlation 
functions measuring how strong the coupling matrix elements correlate at different collective 
deformation parameters q{t) and q{t') , and Q{E) is the nuclear many body level density at 
excitation energy E. 

II. ENERGY RATE 

In order to define properly dynamics of the classical collective parameters q{t) within the 
cranking approach, one has to consider total energy of the nuclear many body system, which 
can be written as 

m = EM+TT{Hm]m}- (5) 

Differentiating over time both sides of Eq. ([5]), we get 

^ = E^^^+E?^2:^ Pn^ + E^n^ + E^^E^ Pnn. (6) 



nn 



The first term in the right-hand side of [6]) describes a change of the collective potential 
energy. The second term depends on the non-diagonal component of the density matrix 
pnm{t) ■ Its time evolution is caused by the virtual transitions among adiabatic states of the 
nuclear many body Hamiltonian. We believe that such a term is a microscopic source for 
the appearance of the collective kinetic energy: 

(,y \ virt 
^Y.&M^v (7) 

where Bij is a collective mass parameter, 

Bij ^ ^ hi,nmhj,mn^nm\.Pnin Pnn]- (8) 
n,m 

The third term in the right-hand side of[6]) is defined by the real transitions between adiabatic 
many body states thus, defining dissipation of energy associated with the nuclear collective 
motion: 

(^) = E ms^, dEEn{E) 11 rft'exp (-^) y^ml mMt') 



d 

""dE 



(9) 
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see Eq. It can be shown that the fourth term in the right-hand side of [6]) is described by 
the slopes of the adiabatic eigenenergies dEn/dqi and within our random matrix approach 
is neglected compared to the other terms. 

Collecting Eqs. ([7]) and one obtains for ([H]), 



where quantities Fi mean forces acting on the collective subsystem given by 

f^f] r+oo ft ( |7-_+'|\ 



(10) 



dqi 



n/T 



d_ 
BE 



n{E) 



d-p{E,t') 
dE 



III. EQUATIONS OF MOTION FOR THE COLLECTIVE PARAMETERS 



To get equations of motion for unknown col 
contributions to the energy rate ( ITOl) are zero 



ective parameters, we assume that all partial 



E, = 0. 



(12) 



Thus, collective dynamics satisfy the following set of equations 



+ 00 



dEEn{E)Y, / dt'(iW 



\t-t'\ 



h/T 



d_ 
'dE 



n{E) 



dp{E,f) 
dE 



.(13) 



can demonstrate that for the constanttemperature level density, Q{) = const ■ exp{E/T), 
where T is the temperature of the nucleus, we obtain a closed set of equations of motion for 
the collective parameters: 

\t-f\' 



f)f] Q.. /■ + 00 



dqi 



T Je. 



dEEVt{E)Y, / c/t'exp 



h/T 



Mm^^it')q,{t') 



X 



_d_ 

dE 



n{E) 



dp{E,t') 
dE 



.(14) 



It should be pointed out that Eq. (fn|) describes the ensemble averaged collective dynamics, 
i. e., averaged over many different random realizations of the intrinsic nucleonic subsystem. 
In this way, of course, we loose information about the quantum fluctuations of the intrinsic 



degrees of freedom of the nuclear many body system which, in principle, may be impor- 
tant. In order to try to take them into account somehow, we introduce a phenomenological 
stochastic force term, ^i(t), into the collective equations of motion (IT^ and requiring that 
the fluctuation-dissipation theorem is hold, 

this, the collective dynamics gets a form of the non-Markovian Langevin equations of motion 

dp{E,t' " 



_d_ 

""dE 



n{E)- 



dE 



IV. NUMERICAL CALCULATIONS 



Now we turn to the numerical determination of nuclear fissions characteristics. We study 
the symmetric fission of highly excited heavy nuclei whose space shape may be obtained 
by rotation of some profile function W'^{z) around z-axis. It is considered a 2-parametric 
family of the Lorentz shapes 

w'{z) = {z'~Q{z'+a)n (17) 

where the multiplier Q = — ColCo/^ + Cl) guarantees the conservation of the nuclear volume. 
Here, all quantities of the length dimension are expressed in units of the radius Rq of the 
spherical equalvolume nucleus. The parameter (o defines an elongation of the figure, while 
the parameter (2 is responsible for the neck of the figure. Thus, in the case of C2 = 00 we 
have a spheroidal shape and for 00 < ^2 < the neck appears. 

The adiabatic collective potential energy of deformation Egg were taken from j^. The 
equations of motion ( fT6l) . ( ITSl) were solved numerically with the help of the simplest Euler 
method with the initial conditions corresponding to the saddlepoint deformation and the 
initial kinetic energy E]^i^ = 1 MeV (initial neck velocity (2 = 0). The numerical calculations 
were performed for the symmetric fission of a nucleus ^se-y temperature T = 2 MeV. We 
define the scission line from the condition of the instability of the nuclear shape with respect 
to any variations of the neck radius: 

d'E„ 



dp 



,2 

neck 



0, (18) 



where pneck = (2/^/(0 (Co /5 + Cl) is the neck radius. 

We considered a time of the nuclear descent from the saddle point to the scission ( |T8l) . 
The total number of 2 x 10^ trajectories Co(^)?C2(^); generated by the different random 
realizations of the random forces ^i(t), were taken into consideration in order to define the 
dynamic path of the system f|T6l) . f|T5|) . To simplify the problem and clarify the role of the 
random force in the nonMarkovian dynamics of the system, we stopped running trajectories 
Co{t), C2(t) when they cross the line Pneck(Co, C2) = Pned*) where Pneck* is the neck radius value 
determined from the Newtonian nonMarkovian dynamics ( i. e., when the stochastic terms 
are absent in Eqs. (fT6|) . (ITS!) ). 

In Fig. 1, we show the Langevin (solid line) and Newtonian (dashed line) dynamic 
trajectories of the neck radius Pneck('^)- Horizontal line in the figure gives the scission value 
of the neck radius Pned* derived from the Newtonian calculation. 

A histogram of the distribution p of the saddle-to-scission times tgc is given in Fig. 2. We 
found that the most probable and mean value of the descent time are significantly smaller 
than the Newtonian estimation of tgc shown by a vertical line. In fact, the random force 
speeds up the process of descent from the fission barrier. Indeed, the action of a random 
force, in some sense, "shakes loose" the system giving rise to a smaller time of motion 
between two given points comparing to the corresponding time for the unperturbed of the 
system. It should be pointed out that this is hold both for Markovian Langevin dynamics 
(it may be demonstrated analytically for some quite simple models) and non-Markovian 
Langevin dynamics. The last feature is demonstrated in Fig. 3, where the mean value (tsc) 
of the descent time is plotted versus the memory time r = h/Tij for the Langevin (solid 
line) and Newtonian (dashed line) paths of the system. 

As seen in Fig. 3, the difference between the two nonMarkovian calculations of the mean 
saddle-to-scission time grows with the memory time. This fact may be explained by the 
correlation properties of the random forces in the nonMarkovian Langevin equations, see 
Eq. f|T6|) . As can be seen from Eq. f|T5|) . the increase of the memory time r amplifies the 
correlation in the collective coordinates Co(''^); C2(^) at two subsequent moments of time t and 
t + At. As a consequence of that, one might expect that a quite big change, occurred for 
the coordinates Co(^)5 C2(^) at time t, will give rise to a sufficiently big change for Co(^)) C2(^) 
even at the next time step t + At. This tendency will be stronger as far as the memory time 
r grows up. Therefore, in average the system will reach the scission faster as compared to 
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the non-Markovian motion of the system. 



V. SUMMARY AND CONCLUSIONS 

Within non-Markovian Langeven approach, we have demonstrated a consistent descrip- 
tion of nuclear large amplitude dynamics, including the memory effects and the random force. 
We have averaged the intrinsic nucleonic dynamics ([2]), over suitably chosen statistics of 
the randomly distributed coupling matrix elements and the energy spacings. Owing to this 
procedure, we have derived the diffusion-like equation of motion for the smeared occupancies 
p{E, t) of the adiabatic many body states. Note that the obtained equation of motion for 
p{E, t) is the non-Markovian one, where the memory effects depend on the width of the 
randomly distributed matrix elements. 

Applying the ensemble averaging procedure, we have also derived the collective mass 
parameter and the internal energy rate. Finally, we derive a set of coupled dynamic equa- 
tions for the macroscopic variables which take into consideration the time variation of the 
occupancies of the intrinsic nuclear states. Following the fluctuationdissipation theorem, we 
have incorporated also the relevant random force into the macroscopic equations of motion. 
The main contribution to the rate of dissipation energy is due to the jump probabilities 
leading to a rate of dissipation which depends essentially on the total energy of the nucleus. 
The final result shows that a time irreversible energy exchange between the collective and 
internal degrees of freedom is possible when the level density increases with energy. We 
have applied our approach to the description of the descant of the nucleus from the fission 
barrier. We show that the random force accelerate significantly the process of descent from 
the barrier for both the Markovian and nonMarkovian Langevin dynamics. We have ob- 
served that the difference between the two nonMarkovian calculations (see Fig. 3) of the 
mean saddle-to-scission time grows with the memory time. This fact may be explained by 
the correlation properties of the random force in the non-Markovian Langevin equations. 
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Figure captions 

Fig. 1. Typical Langevin (solid line) and Newtonian (dashed line) trajectories of the 
neck radius Pncck of the system flTBl) . f[T^ . Horizontal line is the value of neck radius at the 
scission (1181) pj^^^^* derived from the corresponding Newtonian dynamics. 

Fig. 2. Histogram of a probability density p of the descent times tsc for the non-Markovian 
Langevin dynamics f|T6l) . f|T5|) at the memory time r = ^f^ij = 2 x 10^^^ s, when the size of 
memory effects in the system is quite small. The vertical line gives the Newtonian estimation 
for the descent time. 

Fig. 3. The mean value of the descent time (tsc) for the non-Markovian Langevin dynamics 
f|T6|) . f|T5|) versus the memory time r is shown by solid line. The Newtonian calculation of 
the descent time is given by the dashed line. 
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